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Abstract 

Experimental data obtained for PbFex/2Ta]^/203 (PFT) and NaNbOs: 
Gd (NaNbGd) single crystals show a diffused dielectric permittivity 
peak that is inherent to relaxor ferroelectrics. However some devia- 
tions from the normal relaxor properties were also observed and are 
under discussion. One of the key features of the relaxors is the exis- 
tence of the Burns temperature, at which the polar regions appear. 
We found out that PFT shows this feature but the NaNbGd dielectric 
behavior is different. We analyse these properties within a phenomeno- 
logical theory and propose a microscopic model. 
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1 Introduction 



Relaxor properties of heterogeneous ferroelectrics are still intriguing although 
the main phenomenon, the diffuseness of the temperature dielectric-permittivity 
peak, was discovered long ago [||]. It is now established that one of the im- 
portant features of relaxors is the polar-regions formation and these regions 
grow with decreasing temperature and reach their maximal sizes at the fer- 
roelectric or glass- type phase transition p|. 

New information, rather useful for understanding the relaxor properties, 
was recently obtained from neutron scattering which establishes the existence 
of a soft ferroelectric mode above the Burns temperature T^, and below the 
freezing temperature Tg{Td > Tg), 0, |[. In the intermediate temperature 
range, the "waterfall phenomenon" is seen |Q which reflects the absence of 
transverse optical mode below a definite wave vector. Another new idea 
stemmed from the diffuse scattering of neutrons in relaxors suggesting a 
uniform displacement of the ions in polar regions along with optical displace- 
ments 0. We use these findings to establish the order parameter in the 
relaxor state. 

The present study considers the dielectric behavior of PbFei/2Tai/203 
(PFT), and NaNbOaiGd (NaNbGd) crystals which have some common fea- 
tures (the diffuseness of the temperature dielectric peak) as well as some 
differences. In PFT, the extrapolated Curie-Weiss temperature, Tqw, is 
above the temperature of dielectric permittivity maximum, T^, hinting that 
a ferroelectric phase transition would occur if the crystal were homogeneous. 
The heterogeneities in the crystal prevent the appearance of the ferroelec- 
tric phase and lead to changes in the temperature dependence of dielectric 
permittivity characteristic of the relaxor state below Tew We connect 
these changes with the appearance of a new order parameter corresponding 
to a nanodomain structure and show an example of such a structure within a 
model considering a heterogeneous material consisting of slabs (for the sake 
of simplicity) having different elastic and dielectric properties. In contrast 
with PFT, we find that Tc^yis lower than and another order parameter 
appears at in NaNbGd. 
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2 Experimental 



PbFei/2Tai/203(PFT), and NaNbOaiGd (NaNbGd) crystals were grown by 
the flux method. The crystals obtained had an isometric (edge dimensions 
0.5 -2 mm) form with the sides parallel to the {100} planes of the perovskite 
prototype lattice. The details of the crystal preparation and characterization 
have been described elsewhere 0]. 

The dielectric studies were carried out in the 10^ -lO^Hz range in the 
course of both heating and cooling at the rate about 2 -3 K/min with the 
aid of the R5083 and E7-12 capacitance bridges. Aquadag electrodes for 
dielectric measurements were deposited on the opposite faces of the as-grown 
crystals. 



3 PbFei/2Tai/20, 



At high temperatures PFT has the cubic perovskite structure in which the 
B-perovskite-position is occupied by the approximately randomly distributed 
Fe^^ and Ta^"*" ions. This random distribution influences the displacements 
of the Pb ions and can result in the appearance of polar clusters whose inter- 
actions and growths can trigger a phase transition into a glassy or ordered 
phase. Fig. 1 shows the real part of the dielectric permittivity of the PFT 
crystal measured at different frequencies in the 10^ -lO^Hz range. The strong 
frequency dispersion of e observed at temperatures well above Td is likely to 
be due to conductivity [^. In addition to the known e(T) curves for PFT 
crystals ^ and ceramics [l], ^, besides a diffused and frequency-dependent 
e(T) maximum typical of relaxors, an inflection is observed in the s(T) curve 
at temperature Tjp, which is 30-40 K lower then T^. The frequency dis- 
persion of e is less pronounced at temperatures below this inflection than 
above it. Such behavior is similar to that observed in PbSci/2Tai/203 and 
PbSci/2Nbi/203 |TD[ at the spontaneous transition from the relaxor to nor- 



mal ferroelectric state (see also [[T^]), but for PFT, it can be also connected 
with a glass-type phase transition which appears due to percolation of polar- 



ization [12, 13], or antiferromagnetic phase transition ||T], ^ (recent results on 



Rietveld refinement of the neutron diffraction data |jTj] does not show any 
symmetry changes in PFT from room temperature down to 10 K when the 
superstructure due to the antiferromagnetic ordering is taken into account). 
The fit of the Curie- Weiss law to the high temperature side of the s(T) 
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curve measured at a high frequency (Fig. 2) provides the extrapolated Curie- 
Weiss temperature Tew ~ 310 K which is substantially higher than T^. In 
the same figure we show the difference between the experimental 1/e curve 
and the Curie- Weiss fit. It is seen that this difference has two portions where 
it behaves approximately linearly with temperature. One of the portions in- 
tersects the temperature axis at temperature which approximately equals 
Tqw: while the other, seen at lower temperatures, intersects the T axis at 
a temperature which is close to the value Typ 246 K obtained from the 
Vogel-Fulcher fit of the dependence on frequency. 

The fact that ~ Tew hints to identify a new order parameter appearing 
at this temperature with local polarization. It is consistent with the idea that 
in PFT, at some temperature, T^, there appear polar regions (T^ is higher 
than Tew because of the diffusness of the phase transition that we will discuss 
in detail in Section 5). At the same temperature the dielectric permittivity 
starts deviating from the Curie- Weiss law. 



It was rigorously shown in Ref. [15| that, in the heterogeneous media 



consisting of a ferroelectric slab and dielectric layers (dead layers) the ferro- 
electric region gets broken into 180° domains in order to decrease the depolar- 
ization field (see Section 6 where we consider the rigorous solution for similar 
inhomogeneous problem: two ferroelectric slabs separated by dielectric inter- 
facial layers). On this basis we propose the appearance of such (possibly 
short-range) randomly oriented domain structure in the relaxor state below 
T(i. We will discuss the consistency of this idea with other experiments in 
Section 7. 



4 NaNbOgiGd 

NaNbOs exhibits an antiferroelectric phase transition at 630° K [|l], [16| 



When being doped with Gd at small concentrations the temperature of this 
phase transition decreases and the width of the temperature hysteresis in- 
creases. We discussed these trends in Ref. |TB| assuming the (1-x) NaNbOs- 



(x)Gdi/3Nb03 solid solution formation. At a; ~ 0.12 the thermal hysteresis 
abruptly disappears and a diffused hysteresis-free dielectric peak remains at 



higher concentrations (Fig. 3). We regarded this behavior in Ref. |T6[ to 
the appearance of a concentration phase transition into the relaxor-like state. 
However there are features which are different in comparison with ordinary 
relaxors. 
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First of all we find that that the Tew found from the Curie- Weiss fit of 
the high-temperature side of the dielectric permittivity is much lower than 
Tm (Fig. 4). This makes scenario described in the previous section not 
suitable for the experimental data for NaNbGd and another scenario should 
be developed. In the next Section we will consider a reason for the diffusness 
of the phase transition. We will show that, in spite of the different meaning of 
the order parameter r] for PFT and NaNbGd, the diffuseness of the dielectric 
permittivity peak in these two cases can have similar origin. 

5 Landau-type theory of the diffuseness of 
the phase transitions in PFT and NaNbGd 



In Ref. [|T6| we considered the strain and average square of local polarization 
as concomitant order parameters to an antiferroelectric order parameter. The 
addition of Gd results in the appearance of some distorted polar regions. The 
order parameter which is connected with these regions and responsible for 
the deviation from the Curie- Weiss behavior can have different meaning in 
different materials. For example, it can be the average square of polarization 
or local deformation or, if one thinks about the electronic subsystem, the 
degree of electron localization. In the case of PFT this order parameter can 
correspond to the magnitude of the local polarization in the nanodomain 
structure. In NaNbGd it can be an antiferroelectric order parameter. We 
use for this order parameter the notation r/, we assume that it is scalar or a 
tensor of the second kind, which in the product with gives a scalar, and 
we write Landau Free energy expansion: 



F = F^ + \{a + qu)P^ + \(3P^ + \^P^ + \\P^r]+ 

+\Arf + \Bri^ - kur] + cu^ - EP - ((Tq + ct)u + ... ^ ^ 

where a = a{T — Tqw): ^ = b{T — T^). Here > Tqwi uis strain and ctq is 
internal stress (caused by Gd in the case of NaNbGd). 

The equilibrium solutions for this Free energy are simple: 
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1. P = 0, = 

2. P = 0,V = ^ 

_ -A+ ^A^ + ABjku + XP^) (2) 
25^ 

^ a + + Ar] ± \J {a + qu + Xr])^ — 4/^7 

Here u is the strain found from the equihbrium condition (it is always 
finite because of the stress cxo; the electrostriction contribution to the Hamil- 
tonian quP^ simply shifts the Curie- Weiss temperature and changes the value 
of the nonlinearity constant (3). The first phase is stable at high tempera- 
tures. Then, at lower temperatures, the phase 2 becomes stable, with the 
order parameter t], and, at even lower temperatures, the third phase, where 
the ferrolectric order parameter (here we do not pay much attention to the 
difference between the glass-type and ferrolectric solutions as we are mostly 
interested in the relaxor state) coexists with the order parameter t]. The 
values of rj and P can be found from the solution of equations (11.3). 

In order to find the dielectric permittivity one can write the equilibrium 
condition with respect to P: 



(a + qu)P + PP^ + jP^ + XPri = E (3) 
By taking the derivative of this equation with respect to E one can find 

^ SodE eoa + qu + 3(3P'^ + 5-fP^ + Xr] ^ ' 

For the phase (2) in (^ P =0 and Xr] = Xb (^\/x^ + 6 — x) /2B where 
X = [T — Trj), 6 = 2ku/b. In this case (^ takes the form 



eoa + qu + Xb {y/x"^ + 6 - x) /2B 



(5) 



From experimental data we found ^o-^'? which is simply the difference 
between 1/e and the high temperature Curie- Weiss dependence (l/e)^^^ = 
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Eq (a + qu). A fit of tlie expression 8^X7] = e^Xh — x) jlB to tlie 

experimentally obtained 1/^ — (l/e)^^^ for NaNbGd is shown in Fig. 5 
(Notice that the same quality fit and the same results were obtained from 
the fit of to the corresponding experimental curve). We find that this 
fit obtained with e^a = 3.2861 ■ 10'^ ± 3.78 ■ 10"^ Tew = 185.75 ± 1.84 K, 
eoXb/2B = 3.91 ■ 10"^ ± 3 ■ 10"^ K^^ = 308.4 ± 0.6 K and 5 = 1084.4 ± 54 
is rather good and explains the deviation of the dielectric permittivty 
from the Curie- Weiss law by the appearance of the new order parameter 
below and by the existence of local stresses. We have obtained a similar 
good fit for PFT (Fig. 5) at the values: eoXb/2B = 7.16 ■ 10-*^ ± 1.2 ■ 10"^ 
K-\ = Tew = 295.6 ± 1.0 K and 5 = 1450.6 ± 80 K^. From these 
data the appearance of the quadratic temperature dependencies of e(T) at 
the maximal position becomes clear: it appears due to the expansion of the 
square root in expression (^ with respect to x for small x. Hence we regard 
the diffuseness of the temperature maximum of the dielectric permittivity in 
NaNbGd to local quenched stresses and fields produced by the Gd impurities 
in the matrix of NaNbOs and local polar fields produced by the random 
distribution of Fe and Ta at the B-sites in PFT. 



6 A model for structural inhomogeneity 

We propose a model for structural inhomogeneity based on the recent finding 
by Bratkovskiy and Levanyuk [|15| who studied a short-circuited ferroelectric 
slab with dead layers. Consider short-circuited ferroelectric slabs of the width 
/ (for the sake of simplicity we consider two slabs) divided by dielectric inter- 
facial layers with the width (Iq at the boundaries of the whole system (were 
the system has metalic contacts) and with the width d in the middle (see 
Fig. 6). 

We show that the strain can be excluded from the equations by the or- 
dinary minimization procedure. Let us consider the term describing the 
electrostriction effect together with the corresponding elastic contribution: 

Helast 

(r) = Caf3ai3u1,p/2 — qzzajSUajsPzPz whcrc Ua^ir) is a component of the 
strain and Pz{t) is a component of local polarization (at finite values of 
the local polarization one may consider a piezoelectric effect instead of the 
electrostriction one taking into account that (P, + 5P) ^ ^ 2P,5P). 
Minimizing this Hamiltonian with respect to the strain component one im- 
mediately obtains: Uaff = qzzapPzPzl Capap (at finite local polarizations Uap = 
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dzapPz/ Cai3ai3 wheie d\s a piezoelectric coefficient depending on the direction 
of the local polarization). It implies that the strain can be excluded from 
the total Hamiltonian for problems with zero or constant stress. Then, we 
look for a lowest energy solution with possible spatial alternation of polar- 
ization for ferroelectrics which are improper ferroelastics. Detailed solution 
is worked out in Appendix A. 

The macroscopic field energy (see Appendix, expression (A8)) is repre- 
sented in the form 



/m = Fm/LS = Jmo + lAal - E^a^ (6) 

where cxo is the average charge on the ferroelectric-dielectric boundary; the 
coefficients are given in Appendix A. The effective field E* differs from the 
average field E = U/stL by the factor 2e^f/e2zL where U is applied voltage, 
Et is the total dielectric permittivity obtained treating the system as a series 
of the capacitors, is the dielectric permittivity of the ferroelectric slab. 
This factor is small only when /(thickness of the ferroelectric clab) is very 
small. This fact implies that the effective acting field E* conjugated with ctq 
for the thick ferroelectric layers and thin dielectric interfacial layers is largely 
enhanced with respect to the average field that makes it possible to change 
the polarization by very small dc fields. 

On the basis of these results we predict the appearance of the alternating 
domains in the ferroelectric slabs separated by dielectric layers in lamellar 
structures. Perhaps some evidence for this was found in the KTaOa/KNbOa 
superstructure JI^ although this has not been understood well yet. 



The equilibrium value of the macroscopic charge at the boundaries ctq 
corresponds to the minimum of the macroscopic Free energy: (Tq = Ueg/STiLg 
[Sg is the dielectric permittivity of the capacitors corresponding only to the 
dielectric interfacial layers, Lg is the total width of the dielectric layers along 
thefield). Correspondingly, the derivative of this charge with respect to field 
gives 



dao Lsg 



(7) 

dE %-KLg ^ ' 



It is seen that this derivative is large if the total relative length of the 
dielectric interfacial layers LgjL is small. This result is consistent with that 
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obtained in [IE for one ferroelectric slab having dead layers. There can be 
some variations of the dielectric region's widths but only their total relative 
width is important. We point out that this contribution to the dielectric 
permittivity arises only from the domain walls movements due to the voltage 
[/which results in change in total polarization, hence on ctq. 

We now relate the stripe-like solution to experimentally observable quan- 
tities. For a given solution of local polarization, we obtain local strains 
through the minimization of energy with respect to local strain |Q as de- 
scribed at the beginning of this section. The local strains are then used to 
obtain the acoustic-like (collective motions of atoms in the unit cell) atomic 

displacements da (r) through the relation: Uai3 = | + §^ j . For sim- 
plicity, we restrict here to analysis in 2 dimensions. 

In addition to the solution with polarization in the domain perpendicular 
to the plane of the strip (the transverse case T) described in this section, 
we consider a longitudinal (L) case where the polarization is in the plane 
of stripe and along the direction perpendicular to the domain wall. Note 
that the polarization in the adjacent domains is opposite in sign, giving 
180°-type domain walls separating them. The polarization is along (1-1) and 
(-11) directions for the transverse case and along (|^) and (-1-1) directions 
in the longitudinal case, and the strip itself runs along (0) direction. In 
Figure 7, we show contourplots of atomic displacements, local strains and 
local polarization for the two cases. It is clear that the atomic displacements 
in the transverse case are almost constant at the boundaries between polar 
domains and the embedding dielectric. In contrast, we find large accumulated 
local strains at the domain walls for the longitudinal case, making it costly 
in energy. 

We use atomic displacements (ia (acoustic modes) and polarization (optical 
modes) in Eq. (2) of [|T^ to obtain inelastic scattering intensity plots, dis- 
played in Fig. 8. The scattering plots for both T and L resemble those shown 
in Figure 4 of Ref. 0. The atomic displacements corresponding to the 
acoustic modes (inhomogeneous strain) in our analysis corresponds to the 
uniform shifts that Ref. p uses to model their data. Our results that the 
patterns are rather similar for both T and L cases implies that more than 
one kind of atomic displacements or structures could be used to model ex- 
perimental results and that the stripe-like solutions for the inhomogeneity is 
consistent with the experimental observation. 
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7 Discussion & Summary 



We have presented new experimental data on perovskite crystals which show 
diffuse temperature dielectric anomalies. We presented a Landau-type theory 
of this diffuseness by considering local stresses produced by inhomogeneities. 
In the case of PFT the principal finding is that the extrapolated Curie tem- 
perature Tew is larger than Tm and coincides with the temperature T^. Fol- 
lowing the recently proposed idea describing alternating domains in a hetero- 
geneous media consisted of a ferroelectric thin film having thin dead layers or 
an epitaxial film grown on a substrate [|1^, [1^ we propose that similar (but, 
most likely, short-range and random) structures appear in relaxors below T^. 
We provided a rigorous solution for an extended heterogeneous media, two 
ferroelectric slabs divided by dielectric interfacial layers, and obtained that 
this lamellar structure also must have alternating domains in the ferroelectric 
slabs. We showed that the contribution of the domain wall movements into 
the dielectric permittivity in such a system is governed by the relative ratio 
of the total width of the dielectric layers and does not depend on the distribu- 
tion of the widths over the layers. We also computed strains, displacements 
and inelastic scattering intensity of such a domain structure which are in line 
with recent experiments ^ ^, The assumption that there are only 
two main kinds of the regions, ferroelectric and dielectric, is consistent with 
the experimental finding showing the existence of two different transverse op- 
tical modes corresponding to a two-band structure of heterogeneous binary 
compounds [§]. 

The alternating domain structure proposed is consistent with many ex- 



periments. Indeed, for instance, the "waterfall phenomenon" 0, could 
be related in this case to (nanoscale) domains, the domain size determines 
the critical wave vector below which the waterfall phenomenon is observed. 
The coupling between local strain and polarization in alternating domains 
leads to atomic displacements consistent with the inelastic scattering results 
and the observation of strong interaction between acoustic and optical 
modes |20[. A strong coupling between the ferroelectric fluctuations and 
domain wall fluctuations leads to strong damping of ferroelectric fluctuations 
below the critical wave vector. 

We have shown that in the model of the domain structure considered 
the contribution of the domain walls movements into dielectric response is 
especially large in the case when the relative width of the dielectric layers in 
the direction of the field is small and it is consistent with the recent theoretical 
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result obtained for a ferroelectric thin film having dead layers [|T^. The 
universal relaxation recently observed in PMN even at temperatures larger 



than Tm ||2T[ can now be understood on the basis of the irreversible domain 



wall movements ||22| (the movements of the domain walls until there exists 
an electric field without further restoration after switching the field off). The 
Debye-type relaxation of domain walls is expected in the cases when there 
is a strong restoring field and this is also consistent with the existence of a 
strong dispersion at Tm [|12], ^ The dielectric permittivity can now be 
found on the basis of the theory considering the phonons coupled to relaxators 
l^, pB, . This coupling in addition to the domains' freezing can explain 



the experimentally evidenced deviation from the Arrhenius law [25, 26, 27 



Some contribution to the dielectric permittivity from the Maxwell- Wagner 
mechanism is also possible. 

Formation of the proposed alternating domain structure is not hampered 
by the depolarization field which is unfavorable for the formation of lone 



polar regions in a dielectric media [28|. The cooperative appearance of the 



alternating polar regions make them energetically more favorable than the 
lone separated polar regions [0. We should stress that the lone polar re- 
gions (for example in chemical clusters) can also exist and contribute to the 
dielectric permittivity in the presence of free electrons that compensate the 
depolarization fields [^]; the change of polarization in this case is caused 
by both the local ionic dipole reorientations |^ and by the compensating 
charge movements. 

A picture emerging from the data analysis is the following. At Td there 
appear precursors of the alternating nanodomain structure with a short range 
alternation of the local polarization. The maximal possible wave vector of 
this domain structure just corresponds to the critical wave vector seen in 
the "waterfall" phenomenon The dielectric permittivity of relaxors 

is high due to the large ferroelectric fiuctuations and due to domain walls 
movements (see also |2^, |2^, |31|, Q). At some temperature, known as the 



freezing temperature, the nanodomains become very large most probably due 
to the polarization percolation []T2| and are frozen thereafter or a ferroelectric 
phase transition can take place. This scenario seems to be valid in the case 
of PMN and other classic relaxor materials as well. 

In the case of NaNbGd Tew lies below Tm, and the presence of the dif- 
fused temperature dielectric permittivity maximum, in our opinion, is due to 
the local stresses and local fields produced by the Gd impurities. The large 
difference between and Tqw iii this case hints that the new order param- 
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eter appearing at has nothing to do with ferroelectric domain structure, 
also resulting in much smaller permittivity. 
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Appendix A 

The Poisson equation inside each of the slabs can be written in the form 



Here 61,63 are the dielectric permittivities in the dielectric layer at the 
boundary of the system and in the middle respectively; £2z,£2x,£2y are the 
dielectric permittivities in the ferroelectric slab in the corresponding direc- 
tions; ipi is the potential in the i-th slab, and ip'-^ are the first and the 
second derivatives with respect to a. The boundary conditions between the 
slabs at the voltage U are: 



^2.<^L(-//2) - eiy^[,{do) = -47ra; S3y^'.,^{-d/2) - £2.^L(//2) = 47r(T 
£2.<^L(-//2) - £3<^33.(^^/2) = -47ra; £i<^L(^o) - £2.^^41. (//2) = 47ra 
<^i(0) = -U/2; Mdo) = M-ff^y, V2Um = M-d/2y, 
<^3(rf/2) = <^4(-//2); <^4(//2) = <^5K); ipm = U/2 



(9) 



where the general solutions are given by 
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. , ^ sinhfciz ^ cosh.(kiz) ^ 

<Pl{Z: Qx: Qy) = Cn . , , , + Ci2 , . < Z < do 

smnfcicio cosh (ki do) 

sinh(A';22;2) cosh(A;22:2) . ^ ^ r /r, 

Mz., <lx, %) = C,,--^kJj2) + ^^^ cosh(fe2//2) ' < < 

sinh(/c3^3) 

^3(^3, g.) = ^31 -j^^^, -d/2 <Zs< d/2 

9^40^4, = ^41 ■ wl , + ^42 . A < ^4 < //2 

smn(/c2//2) cosn(/c2//2) 
sinh fci^^s cosh(A:i2:5) 

</:'5(^5, = ^51 ,,, + C'52— TTT-TT, < < do 

blllll ft 1 uo CObll yrolOjQ) 

(10) 

Here 2:2 = 2; + do + //2, -23 = Z2 + f /2 + d/2. The interfacial charge 
can be represented as a wave corresponding to alternating domains in the 
X and y directions placed in the checkboard (it is close to the honeycomb 
structure) or stripe-type style: 



2P 

a = (1 - e*«-"-i) (5g^o, 7^ 0, stripes (n) 

AP 

(7 = ' (i_e^g^«^i)(l-e^'?''"i^), g^ + g^^O; checkboard 

Qx-'- xQy-'-y 

Here P5 is the local polarization magnitude inside domains; at zero wave 
vector 0" = (To where (Tq can be found from the equilibrium condition which 
we will derive below. The lengths of the domains in the xand ^directions 
{ttai and aa2 where a = x,y for polarization up and down, = Oqi + 0^2) 
depend on the electric field magnitude. 

Above we considered only the case when the polarization in the ferroelec- 
tric slabs is parallel, as, in the antiparallel case, the energy of the system 
is larger. The solution of the above equations in this case has the following 
symmetry conditions: C42 = -C22; C41 = C21; C52 = -C12; C51 = -Cu- 
Prom the equality of the potentials at the interfaces (see conditions (A2)) one 
obtains: C31 = -C21 - C22; Cn = -C21 + C22 - C12; Cu = at ^ g^ _^ 
and C12 = —U/2 at q'^ + Qy = 0. From the Poisson equations in (A2) we have 
C21F1 — C22F2 — CuGq — —Attu and C31G1 — C2i-^i — C22F2 — Ana, where 
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Go = £1^1 coth (fci^o), Gi = £3/^3 coth (A;3(i/2), Fi = £22^2 coth (A;2//2), 
-^2 = £2z^2 tanh (A;2//2). Excluding C3iand Cu one has:(Fi + Go)C'2i — 
(F2 + G'o)C22 = -47rcr - C12G0, -(Fi + G'i)C2i - (F2 + Gi)C22 = 47ra. Under 
these conditions we obtained the following solution (only one coefficient will 
be necessary for the Free energy calculation) for ql + 7^ 



47ra (2F2 + Go + Gi] 



21 



(Fi + Go)(F2 + Gi) + (Fi + Gi)(F2 + Go) 



(12) 



and for q^ + ql 



0: 



Lg , ^tf 



C„.^-4..-^ + Cj^ (13) 

Here L/e^ = Lg/Sg + L//£22, ^g/e^g = d/e-i + 2c/o/£i, = 2/. 

The fci, A;2 and can be found from the Poisson equations: k\ = = 
ll + lyy ^22^2 = ^^xQI + ^2yQ'y- Finally one obtains the Free energy in the 
form: 



F = Fo-2Sj2 C2Mq., qy) - (14) 
qx,qy 

where 5* is the surface area covered by the domains. The ffist term includes 
the Landau expansion with respect to polarization, domain wall energy and 
the energy connected with the appearance of polarization under the bias field; 
the second term describes the electrostatic energy of the boundaries; and the 
last term describes the capacitor energy due the bias field. 

From (0) the macroscopic field energy (the terms corresponding to + 
— 0) can be represented in the form 



fM = Fm/LS = Jmo + lAa^ - E*a^ (15) 

where /a/o = StU^ML^ A = ML'^ {e^./Lf + eg/Lg)-\ E* = etLfU/e2^L\ 
The alternating domain period can be now found from the equilibrium 
condition for the alternating field contribution to the Free energy 
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NfASP^ ^ SvrSa^ (2F2 + Go + Gi 

J^alt — 1- / 



(Fi + Go) (F2 + Gi) + (Fi + Gi) (F2 + Go) 

(16) 



where =4 in the case of the checkboard order (a = axi = ax2 = o-yi = CLy2) 
and =2 in the case of the striped domains (a = a^i = ax2] CLyi = o-yi = 00); 
A is the characteristic domain wall width. 

Especially simple solution can be obtained in the case when the dielectric 
and ferroelectric layers are thick in comparison with the domain width [15,34] 
and e2x = £2y- 



l NfAie2g + €^){e2g + es) 

64^ (2^29 + £1 + £3) ^ ' 

where e2g = y/£2z£x', ^ does not depend on aand for the checkboard order we 
have: 



4 



00 00 



^0.813 (18) 

i=o j=o (2i + 1)' (2j + 1)' ^J {2i + if + (2j + 1)' 



In the case of the striped domains [|34 

TT C 1 , s 

— = > ^ ^ 1.0518 (19 

2 t^(2z + l)^ ^ ' 

The domain width increases with dielectric permittivity and the thickness 
of the ferroelectric slab. 

In the case of thin dielectric interfacial layers and comparatively thick 
ferroelectric slabs in the stripe-type order and at ei = 63 = ^Je2z^2x■, d/2= 
do the problem can be reduced to the case considered in Ref. JTSf, and the 
final result is 



a = 0.95rfexp ( 0.4^ ) (20) 
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where qk is the Kittel domain size parameter which coincides with (P!7[ ) . This 
expression shows that the domains extremely rapidly grow when the relative 
size of the dielectric interfacial layers in the direction of field becomes thinner. 
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Figure captions 

Figure 1. s'(T) dependencies for the PFT crystal, measured at different 
frequencies. 

Figure 2. Temperature dependencies of e' (1) and 1/e' (2) measured 
at lO^Hz for the PFT crystal. Curve 3 shows the difference between the 
experimental 1/e dependence and the Curie- Weiss fit oil/e . Straight solid 
lines are guides to the eye. 

Figure 3. Temperature dependencies of e' measured at 10^, 10^^,10^, 
lO^Hz. for 0.88NaNbO3-0.12Gdi/3NbO3 crystal. 

Figure 4. Temperature dependencies of e and 1/e measured at lO^Hz for 
0.88NaNbO3-0.12Gdi/3NbO3 crystal and the difference between the experi- 
mental 1/e dependence and the Curie- Weiss fit of 1/e. Straight solid lines 
are guides to the eye. 

Figure 5. The fit of the model expression to the experimental data. 

Figure 6. The model lamellar structure consisted of ferroelectric and 
dielectric slabs. 

Figure 7.Contourplots of (a) atomic displacements dxir), dy{r), (b) local 
strains Uxxij), Uyy{r), (c) Uxyij) and local polarization Pxir) for the alter- 
nating domain-stripe solutions of the model inhomogeneous ferroelectric in 
2-dimension. 
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Figure 8. Contourplot of the inelastic scattering intensity for the alter- 
nating domain-stripe solution discussed in Section 6. 
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